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We generalize the oscillator model of a particle interacting with a thermal reservoir 
by introducing arbitrary nonlinear couplings in the particle coordinates. The equilibrium 
positions of the heat bath oscillators are promoted to space-time functions, which are 
shown to represent a modulation of the internal noise by the external forces. The model 
thus provides a description of classical and quantum dissipation in non homogeneous envi- 
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multiplicative noise and a position-dependent fluctuation-dissipation theorem associated to 
non homogeneous dissipative forces. When time- modulation of the noise is present, a new 
force term is predicted besides the dissipative and random ones. The model is quantized 
to obtain the non homogenous influence functional and master equation for the reduced 
density matrix of the Brownian particle. The quantum evolution equations reproduce the 
correct Langevin dynamics in the semi classical limit. The consequences for the issues of 
decoherence and localization are discussed. 
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1. Introduction and summary. 

In the Langevin approach to classical Brownian motion |lj the total system, globally 
isolated, is divided into two parts, the central particle and the environment. The effects of 
the environment on the otherwise free particle are taken into account by adding two force 
terms in Newton's equation, a deterministic damping force —M^ x and a random force 
R(t) (for simplicity we consider only one-dimensional motions, since the generalization to 
higher dimensions is straightforward): 

Mx = -M^ x + R{t). (1) 

The former expression is the Langevin equation, where M stands for the mass of the 
Brownian particle, 7 is the damping constant, and R(t) is a stationary Gaussian Markov 
process. The expectation value (R{t)) of the random force is assumed to be zero, 

(R(t)) = 0, (2) 

so that the macroscopic equation for the expectation (x(t)) coincides with the classical 
equation of motion of a damped particle, M'x c i = —M^ x c [. The additional requirement 
that, for t — > oo, the particle reaches thermal equilibrium with the environment, so that 
M(x 2 )/2 = k B T/2, requires that: 

(R(t)R(s)) = 2 lo Mk B T5(t - s), (3) 

where k B is Boltzmann's constant and T is the absolute temperature. R(t) is thus a white 
noise. Eq.(3) expresses the well known fluctuation-dissipation theorem for the particular 
case of white noise, in which the correlation function Kit — s) oc 5(t — s). 

However, there is no Lagrangian and therefore no action principle allowing to derive 
the phenomenological equation of classical Brownian motion, eq.(l). This shortcoming 
of the classical theory makes it very difficult to build a quantum theory of Brownian 
motion and to treat the quantum mechanics of open systems because canonical or path- 
integral quantization require the knowledge of the microscopic dynamics of the system, 
either Lagrangian or Hamiltonian. 

A possible way to overcome this difficulty is to introduce a suitable mechanical model 
of dissipation. A particular one, the so-called oscillator model, has grown to a privileged 
position since its inception in the early sixties (f2|, |§, (§]). The oscillator model, in its 
classical version, provides an appealing alternative derivation of the Langevin equation of 
classical Brownian motion. Upon quantization, the model provides a theory of quantum 
Brownian motion which, in the semiclassical limit, reduces to the known phenomenological 
Langevin dynamics. 

The first step towards a mechanical model of homogeneous classical and quantum dis- 
sipation is to consider the central particle and the dissipative environment as forming a 
large isolated dynamical system. 

The second step is to consider the dissipative environment as composed by an infinite 
collection of independent harmonic oscillators {q n }, with frequencies {u n /2ir}, linearly 
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coupled to the "central" particle x of mass M. Then, if suitable assumptions are made on 
the statistical distribution of the initial conditions of the oscillators, the central particle 
exhibits Brownian motion with homogeneous noise in the limit of continously distributed 
frequencies. 

The Lagrangian of this many-body mechanical model is 

L{x,x,{q n ,q n \,t) = —x + ^ \^q n ~ -^u n {q n - x) j , (4) 

where for simplicity all the oscillators are assumed to be of equal mass m. The first term 
in the r.h.s. of eq.(4) is the free Lagrangian Lq of the central particle, while the remaining 
part will be referred to as the oscillator Lagrangian L osc ; it describes the whole environment 
and its interaction with the central particle. The oscillator Lagrangian actually describes 
a field-particle coupling, for any free field can be essentially viewed as an infinite collection 
of harmonic oscillators. 

The oscillators are assumed to be in thermal equilibrium at the initial time t according 
to a canonical distribution: 



(<?io - xo) 


= o, 


(fto) 


= o, 


muJiUjUqio - x )(q j0 - x )) 


= k B TSij 


m(qi q j0 ) 


= k B T8ij 


(<?io(<?jo - ^o)) 


= o, 



(5) 



where dij is the Kronecker symbol and qio = qi(to), <jio = (ft (to), and xq = x(t ). The 
classical homogeneous Brownian dynamics of the central particle is obtained by eliminating 
the oscillators coordinates from the Euler-Lagrange equations. The resulting equation of 
motion for the central particle x is the generalized Langevin equation with homogeneous 
noise, considered by Mori and Kubo in the context of nonequilibrium statistical mechanics 

Mx = - [ dsK(t-s)x{s)+R(t). (6) 



The homogeneous noise R(t) and the homogeneous memory kernel K(r) are, respec- 
tively, 

R(t) = mu l \ (ino ~ x ) cos[^ n (t - t )] + — sin[u n (t - t )} \ , (7) 

n L LO n ) 

K i T ) = mul cos (u} n r). (8) 
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Note that the initial oscillator coordinates and velocities are just the Fourier coefficients of 
R(t). By using initial conditions (5), it is straightforward to show that R(t) is connected 
to the memory kernel K(t) by the fluctuation-dissipation theorem: 

(R(t)R(s)) = k B TK(t - s). (9) 

The form of the correlation function depends on the oscillator distribution. For instance, 
white noise corresponds to an oscillator distribution function G(uj) = 2M r y Q /'Kmuj 2 . In this 
case the correlation function is 

POO 

K(t)= duG(u)muj 2 cos (out) = 2 7o 5(t) (10) 
Jo 

and eq.(6) and (9) reduce to eq.(l) and (3) respectively. 

It should be remarked that the dissipation-fluctuation theorem was originally obtained 
on grounds of internal consistency of the Langevin equation ||, while in the oscillator 
model it is a natural consequence of the form of the Lagrangian (4) and of the initial 
conditions (5). 

The model Lagrangian (4) is the only quadratic one which is both reflection and trans- 
lation invariant. These invariance properties must be shared by any model Lagrangian 
describing a globally isolated system, in order for eq.(l) to be invariant. 

In fact, the homogeneous oscillator model can be effectively visualized as a set of masses 
m bound to the mass M of the central particle by linear forces ||. In particular, the 
coordinate x of the central particle represents the equilibrium position of the oscillators. 
This explains the presence of xq in the initial data of the oscillators, eq.(5). Lagrangian 
(4) has been first proposed in ref. ||, and it has been later considered in different contexts 

(0,1,0,0)- 

Models with an oscillator Lagrangian which is not translation and reflection invari- 
ant produce non physical infinite force terms and ill-defined quantities which depend on 
frequency cutoffs, both in the classical and in the quantum regime (P,|ll|, [|l^] , fl3| ) . 



Consider now the case in which an external potential V(x, t) is present. The usual 
generalization of eq.(6) is 



Mx + d x V(x, t) = - f dsKit - s)x(s) + R(t), (11) 

Jta 



obtained by simply adding the external force —d x V(x,t) = dV(x,t)/dx in the generalized 
Langevin equation. It can be derived from the oscillator model by simply adding the 
potential term V(x,t) in Lagrangian (4). In many physical situations this equation or its 
white noise limit provide a satisfactory description of classical Brownian dynamics. 

However, from a general point of view some inconsistencies arise in this approach. 

At a phenomenological level one can see, by inspection of eq.(ll), that the addition 
of an external potential V(x,t) has affected the mechanical part of eq.(6), on the left 
hand side, leaving unchanged the environment-induced force terms, on the right hand side. 
Instead, a corresponding influence should be expected to appear also on these terms, since 
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the environment surrounding the Brownian particle is made up of particles as well. Thus, 
eq.(ll) cannot be in general a correct description of Brownian dynamics in the presence of 
external fields. 

There are many physical situations in which the non homogeneous nature of the en- 
vironment has observable consequences on the dynamics of the central degree of freedom. 



Some examples are provided by systems escaping an oscillating barrier stochastic 
resonance or stochastic parametric oscillators ||15| . 

The purpose of the present paper is to provide a reliable model for the description 
of classical and quantum Brownian motion in non homogeneous environments. This is 
achieved by introducing a suitable generalization of the classical homogeneous oscillator 
model. 

A true mathematical modellization of the bath and the central particle as forming a 
closed system is really sound as long as V(x, t) = 0. In this case L and L osc , as well as the 
total Lagrangian, share the same invariance properties. If we allow for external potentials 
V(x, t) acting on the central particle, we are in fact considering situations in which the 
whole system (bath + particle) is open and it is then natural to consider it coupled to 
the rest of the universe through external forces which act both on the particle and on the 
thermal oscillators. In this way one treats the central particle and the oscillators on equal 
dynamical footing. 

This approach enables one to study those physical situations where strong enough 
couplings with external fields do not leave the environment unaffected. 

The generalized Lagrangian is introduced in Sec. 2. Besides adding an external potential 
V(x,t), we let external forces F n (x,t) act upon each oscillator q n with coupling strength 
c n , in order to describe the action of the external potential on the environment. This is 
shown to be equivalent to shifting the equilibrium positions of the oscillators from x, the 
position of the central particle, to the positions Q n (x, t) = c n F n (x, t)/mu}^. It is also shown 
that the functions Q n represent a space-time modulation of noise and dissipation. 

We treat the classical model first, and for it we derive the Euler-Lagrange equation of 
motion for the central particle. It turns out to be a modified generalized Langevin equation 
with a colored and multiplicative noise term R(x, t) that can be expressed as the gradient 
of a stochastic potential V(x,t). Further, a position-dependent memory kernel K replaces 
the homogeneous kernel K(r) in the dissipative force term. A new deterministic force 
is predicted besides the random and dissipative ones when time-modulation of the noise 
is present. A generalized position-dependent fluctuation-dissipation theorem is derived 
connecting the memory kernel K and the noise R. The corresponding modification in the 
Fokker-Planck equation is also discussed. 

Further theoretical and experimental effort is needed to verify the predictions of the 
model. To this end, we study in detail the different physically significant limits of the 
Brownian dynamics in the non homogeneous background. In particular it is shown that 
in the white noise limit the friction coefficient in eq.(l) acquires a space-time modulation: 
7 — ► 7(2, t), and that a precise relation between j(x, t) and the stochastic potential V(x, t) 
holds. 

In Sec. 3 the model is quantized via Feynman path integration. The non homogeneous 
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influence functional and the master equation for the reduced density matrix of the central 
particle are obtained in closed form. We discuss the semiclassical limit and verify its con- 
sistency with the results obtained from the classical model. The expressions obtained are 
compared with those already derived in the literature in order to analyze the consequences 
for the issues of decoherence and localization. 



2. The classical model. 

2. a Langevin equation and fluctuation-dissipation theorem. 

Our starting point is the total classical Lagrangian, eq.(4), of a central particle x in 
interaction with an infinite set of harmonic oscillators {q n }. We now consider that a non 
zero external potential V(x,t) acting on the central particle is present. In this case the 
interpretation of L osc as the Lagrangian of the whole environment does not hold and there 
is no more reason to request its reflection and translation invariance. In order to take 
into account the influence of the external potential on the environment, we let the generic 
oscillator q n interact with the external field through a coupling function F n (x(t),t). The 
new total Lagrangian is then 

L(x,x,{q n ,q n },t) = —x 2 - V(x, t) + ^ | — (q 2 n - u 2 n q 2 n ) + c n q n F n (x, t)j , (12) 

where the C n S £1X6 constants. This Lagrangian, already considered in [fTq| , cannot lead to 
the correct generalized Langevin equation. The reason is that in the homogeneous limit, 
that is by choosing the simple linear particle-oscillator coupling F n = muj 2 x/c n , it does not 
reduce to Lagrangian (4), but to the original non invariant Feynman- Vernon Lagrangian, 
where the renormalization potential ^ n mu) 2 n x n /2 is missing. 

In fact, it can be shown that the generalized Langevin equation obtained from the 
Lagrangian of eq.(12) contains a force term which is the gradient of the nonlinear renor- 
malization potential AV = — J2n c nF 2 /2mio 2 . 

In order to avoid the presence of this ill-defined quantity in the equations of motion we 
redefine the Lagrangian by adding the renormalization potential — AV. It is convenient to 
introduce the quantities Q n (x(t),t) = c n F n (x(t),t)/muj 2 , with the physical dimensions of 
a coordinate. Then Lagrangian (12) becomes 

L(x, x, {q n , q n }, t) = yi 2 - V(x, t) + ^ {\ m ol ~ \ muJ l [in - Qn(x, t)] 2 | . (13) 

It is to be noted that formally this Lagrangian and the homogeneous one, eq.(4), differ 
only in that the external potential V(x, t) is present and the oscillator equilibrium positions 
are now represented by the functions Q n (x,t). 
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Lagrangian (13) defines the most general oscillator model which is linear in the oscillator 
coordinates, thus preserving the exact solvability of the oscillator sector, and nonlinear in 
the central particle coordinates. This Lagrangian has been considered by Zwanzig |7| and 



by Lindenberg and Seshadri [l7||, in the classical regime for the case in which the Q„'s do 
not depend on time. 

The equation of motion for the central particle can be derived from Lagrangian (13) 
following the same procedure applied in the homogeneous case. Introducing the compact 
notations x(t) = x, x(s) = y, and t — s = r, one obtains the following generalized Langevin 
equation 



M'x + d x V{x,t) = - I dsK(x,y,r)y - f ds®(x,y,r) + R(x,t) . (14) 

J to J to 



The first and the third terms on the right hand side of the equation are the dissipation 
and fluctuation forces respectively. The second term in the r.h.s. is an additional force, 
neither random nor dissipative. Its physical origin is explained below. 

The expressions for the fluctuating force R(x, t), the memory kernel K(x, y, r), and the 
deterministic kernel $>(x, y, r) are 



R{x,t) = J2 muJ l d xQn{x,t)Uq n0 - Q nQ ) cos[u n (t - t )] + —sm[u n {t - t )}\, (15) 



K(x,y,rj =J2{dxQn{x,t)dyQ n (y,s)}mujlcos(uj n T), (16) 

n 

®(x,y,T) =^2^d x Q n (x,t)d s Q n (y,s)jmulco8(uj n T), (17) 

n 

where Q n0 = Q n {x(t ),t ). 

The oscillators {q n } of the environment are assumed to be in thermal equilibrium 
at time to around their effective equilibrium positions Q n (x,t) according to a canonical 
distribution, in complete analogy with the homogeneous case: 

(to - Qio) = , 

(to) = , 

muiUj((qi - Qio)(qj - Qj )) = k B T5ij , (18) 

m(q i0 q j0 ) = k B Tbij , 

(tofeo - Qjo)) = 0. 
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Comparing the definition (15) of R with the above statistical conditions on the oscil- 
lators, one can prove that the fluctuating force is a Gaussian stochastic process with zero 
mean 

(R(x,t)) = 0, (19) 

and covariance 

(R(x, t) ■ R(y, s)) = k B TK(x, y, r) . (20) 

It is remarkable that a fluctuation-dissipation theorem holds even in this non homoge- 
neous generalization of the oscillator model. Furthermore, no restriction has been assumed 
on the different time scales of the noise and of the external potential. From the definition 
(15) and the statistical properties (19)-(20), it follows that the fluctuating force R is a col- 
ored multiplicative process, in general not factorizable in the product of a time-dependent 
and of a space-dependent part. 

If the functions Q n have no explicit time dependence, then from eq.(17) it follows that 
the last force term in the Langevin equation is absent. On the other hand, if the functions 
Q n do explicitly depend on time but not on x(t), then from eq.(16) it follows that the 
dissipative force is absent, and we are left with a purely time-dependent force, a case 
already considered by Feynman and Vernon . Finally, the case of homogeneous noise is 
regained as V = and Q n (x,t) = x, Vn, and eqs.(13)-(16) and (18)-(20) go over into the 
corresponding eqs.(4),(10),(7),(8) and (5), (2), (9). 

The Langevin equation obtained above is very general because, as we earlier remarked, 
the noise R and the correlation function K do not factorize. However, when all the functions 
Q n are equal to the same function Q(x,t), the equation simplifies into: 

Mx + d x V(x, t) = d x Q(x, t) jf2(t) - J* K(t) (d s Q(y, s) + yd y Q(y, s))ds^ , (21) 

where R(t) and K{r) are those of the homogeneous case, given by eqs.(7)-(8), and, as 
before, y stands for x(s). As a consequence, the generalized fluctuation-dissipation theorem 
(20) reduces to the homogeneous form (9). 

Inspection of eq.(21) shows that the effect of the new oscillator equilibrium position 
amounts to a space-time modulation of the fluctuating force, which can be rewritten as 

R(x,t) = d x (R(t)Q(x,t)) = -d x V(x,t) . (22) 

We are thus naturally led to the interpretation of —Q(x,t) as a deterministic modulation 
of the fluctuating potential V. It is now clear that the last force term on the right hand 
side of eq.(21) derives from the time modulation alone of the stochastic force. Namely, 
the model predicts that if there is a time modulation of the noise, then, besides the usual 
dissipative and random terms, a new (non dissipative) force appears. 

If the functions Q n {x, t) are different from each other, then the random force R(x, t) in 
eq.(14) can still be written as the gradient of a stochastic potential V(x,t), given by the 
expression 

V(x, t) = - muJ n Q n (x, t) I (q n0 - Qno) cos[u n (t - t )\ + — sm[u n (t - t )\ \ . (23) 
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In this case we preserve the same interpretations of the functions Q n . The above 
relation can be seen as a particular Fourier expansion of the random potential, in which 
every Fourier component is modulated in space and time by the functions Q n (x,t). 

In the white noise limit K(t) — > 2M r y 5(r) and eq.(21) reduces to 

Mx + d x V(x, t) = -M7(x, t)x + R(t)d x Q(x, t) - M^ d x Q(x : t)d t Q(x, t) , (24) 

where we have defined the effective friction coefficient as 

7(x,f) =-y (d x Q(x,t)) 2 . (25) 

We see that the nonlinear coupling to the linear bath modifies the white-noise Langevin 
equation by introducing a nonlinear dissipative force, a multiplicative noise term, and a 
deterministic force due to the time-variations of the coupling. The friction coefficient 7 is 
modulated by the square of the function that modulates the white noise R(t), as expected 
from the fluctuation-dissipation theorem (20). For Q n (x, t) = Q(x, t), the following relation 
holds: 

R(x,t) = -d x V(x,t) = R{t)d x Q{x,t) = J^^-R(t). (26) 

V 7o 

It is instructive to analyze the overdamped regime as well. In this case however, one 
must be aware of the fact that the very existence of the overdamped limit rests on the 
hypothesis that the scale of the time- variations of Q(x,t) is slower than the relaxation 
time l/7o- 

Assuming that this constraint is satisfied, in the overdamped regime one has that the 
function modulating the fluctuation goes into its inverse: 

i= -^(j^kQ-^k, +d,Q )- (2?) 

When d t Q = 0, according to the above expression the particle moves under the influence 
of an effective force —d x V/d x Q 2 . Compared to the homogeneous case, the modulation of 
the noise can strongly modify the intensity of the external force in the overdamped limit 
but cannot change its direction. However, when dtQ ^ 0, also the sign of the effective force 
can be reversed if the scale of time- variations of Q(x,t) is small enough. 

Finally, if Q does not explicitly depend on time, eq.(27) suggests a simple method 
to determine d x Q and so the nonlinear coupling Q as well. From a measurement of the 
velocity x and the knowledge of the external potential V(x,t), it is possible to infer the 
non homogeneous nature of the environment and to determine d x Q given that, on average, 
d x Q 2 = -d x V/M lo x. 

2b. Fokker-Planck equation. 

In the case of white noise the Langevin equation 

Mx(t) + d x V(x, t) = -M lo x(t) + R(t) , (28) 



9 



is equivalent to the Fokker-Planck equation 

d t W = -pM~ l d x W + d p (wd x V + KpW) + M^- x d 2 v W 



(29) 



for the phase-space probability density W(x,p, t), where p is the momentum of the particle 
and = ksT. 

In the overdamped regime the above equation reduces to the diffusion equation for the 
configurational probability density P(x,t) = J dpW(x,p,t): 



d t p = (m 7() ) '[^(p^f) + (3~ 1 d 2 x p 



(30) 



Considering now the case of non homogeneous noise, the Fokker-Planck equation equiv- 
alent to the white noise Langevin equation (24) is 

d t W = -pM- l d x W + d p [( 1P + M lo d x Qd t Q + d x v)w] + M 1 (5~ 1 d 2 p W . (31) 

It is to be noted that since the same function j(x,t) appears both in the drift and 
in the diffusion term, the equilibrium (stationary) solution of eq.(31), if Q and V do not 
depend on time and if such a solution exists, is just the Maxwell-Boltzmann distribution 
W eq .(x,p) oc exp [-p (p 2 /2 + V(x))}. 

In the overdamped regime, one obtains the diffusion equation for the configurational 
probability density P(x,t): 

d lQ \ 



d t P = d x 



'd x V p- 1 d t Q 
M 7 M lo d x Q (d x QY 



+ ^ r (d x [(d x Q)- 2 ])d x P. (32) 



M 7o 



This equation can be used to study specific examples of nonlinear couplings with given 
external potentials. It is easy to see that the modifications induced by the non homogeneous 
terms may lead to interesting non perturbative effects. We shall present in a forthcoming 
paper a detailed treatment of the non homogeneous Brownian dynamics in bistable and in 
periodic potentials. 



3. The quantum model 

3. a Influence functional. 

We shall now discuss the quantum version of the generalized oscillator model, by quan- 
tizing the Lagrangian (13) via Feynman path integration, and derive (see Sec.3.b) the 
relevant evolution equations for the reduced density matrix p of the central Brownian 
particle. 

We begin by considering the action S[x, q] of a system formed by a central particle x 
and just one oscillator q of natural frequency u/2ti: 



S[x,q] = S [x] + S osc [q,x] 



£ dt [^X 2 - V(x,t)} + £ dt^ {q 2 - U 2 {q - Q(x,t)) 2 } , 



(33) 
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where t a and t b are two arbitrary initial and final times, respectively. The functions So and 
S osc are respectively the actions of the isolated particle and of the oscillator, while Q(x,t) 
is an arbitrary space-time function. 

For the sake of clarity, we shall assume here and in the following that the initial density 
matrix of the total system factorizes in the product of the density matrices of the two 
subsystems. In the coordinate representation, denoting by /j and j[ two different values 
of an observable quantity f(x,t) at time ti, the former statement reads 

p(x a , x' a , q a , q' a , t a ) = p(x a , x' a , t a ) ■ p*(q a - Q a , q' a - Q a ). (34) 

Some remarks should be made about the form of this density matrix. Equation (34) 
describes the result of a measurement in which the state of the central particle has not 
been affected by the environment (0, [|K|, Hl8|). The density matrix p(x a , x' a , t a ) represents 
the initial state of the central particle, and p*(q a — Qa) q' a ~ Qa) is the equilibrium density 
matrix of a harmonic oscillator at an inverse temperature (3 = 1/kgT with its equilibrium 
position at q a = Q a . It has the expression 



P*(qa-Q a ,q a -Qa) = F((3)exp 



2h 

2(q a -Q a )(q' a -Qa 



((q a -Qa) 2 + (q a -Qa) 2 )coth((3huj) 



(35) 

where F(f3) = y 'mcu / '-kTi coth.{j3huj / '2) is a normalization factor such that / dqp*(q,q) = 1. 

The oscillator density matrix only depends on q a and q' a through the differences q a — Q a 
and q' a — Q a , where Q a , in analogy with the case of homogeneous dissipation |1§, can be 
shown to be given by 

/=, Q( x a, ta) + Q( x a i ^a) loa \ 
Qa = ^ • ( 36 ) 

This expression accounts for the fact that, compared to the homogeneous case, in which 
Qa — {%a J r x ' a )l^'i ^ e equilibrium position of the oscillator at time t has been shifted from 
x(t) to Q(x(t),t). 

The density matrix evolves in time from t = t a to t = through the action of the total 
density matrix propagator J(b\a) = J(x b , x' b , q b , q' b , t b \x a , x' a , q a , q' a , t a ): 

p(x b ,x' b ,q b ,q' b ,t b ) = / dXadx' a dqadq' a J(b\a)p(xa,x' a ,q a ,q' a ,t a ). (37) 



The total propagator J(b\a) is simply given by the product of the two propagators of 

the total wave function: 

r x b r%'h rib ril ( i 1 

J(b\a)= Vx Vx' Vq Vq'exp\-{S[x,q) -S[x',q'))\, (38) 

Jx a Jx' a Jq a Jq' a in J 

where S[x, q] is given by eq.(33). The reduced density matrix of the central particle at 
t — t b is obtained by integrating the total density matrix in the oscillator coordinate, 

p(x b ,x' b ,t b ) = / dq b p(x b ,x' b ,q b ,q' b ,t b ). (39) 
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Following the same procedure as in the homogeneous case (0, fllTl), the time evolution 
law of the reduced density matrix is 

(40) 



(41) 



p{x bl %bi ^b) J dx a dx a J e j f{x\)i Xfo } t b \x a} X a} tg^p^Xa, X a} t a 

where the effective propagator J e ff can be written as 

Jeff(x b ,x' b ,t b \x a ,x' a ,t a ) = / Vx Va/exp[-S[x,x' 

J X a Jx' a \/l 

Here we have introduced the effective action 

S[x, x'] = Sq[x] — Sq[x'] + $[a;, x'), 



(42) 



where $[a;, a;'] denotes the influence phase ||, and must not be confused with the deter- 
ministic kernel <&(x,y,T) defined by eq.(17). Introducing, as before, the compact notation 
x' = x'(t), x = x(t), y' = x'(s), y = x(s), and t — s = r, the expression for the influence 
phase reads 



x'] 



dt \ ds 



+ 



+ 



Q(x', t) - Q(x, t) [d s Q(y, s) + d s Q(y', s) J cos(u;r) 

Q(x', t) - Q(x, t)] {yd y Q(y, s) + y'd y ,Q{y', s)} cos(wr) (43) 

i coth(/5^/2) / dt ds Q{x',t) - Q(x,t) \Q(y' , s) - Q(y, s)\ cos(co>t) . 

2 J t a J t a 



dt / ds 



Because of the dynamical and statistical independence of the oscillators, in the case of 
an environment with infinite degrees of freedom the influence phase can be calculated by 
simply summing up the influence phases of all the oscillators. 

If all the Q n J s are taken to be equal to a certain function Q, the influence phase assumes 
the form 



x'] 



l r tb r* 

dt / ds 



+ 



2 

l r tb 



dt / ds 



+ i dt ds 



Q{x', t) - Q(x, t)] {d s Q(y, s) + d s Q(y\ s) }k(t) 
Q(x', t) - Q(x, t)] {yd y Q(y, s) + y'd y ,Q(y', s)}k(t) 
Q{x', t) - Q(x, t)] {Q(y', s) - Q(y, s)}a(r) , (44) 



where 



a(r) = - moj z n coth(/3fi,u; w /2) cos(cj ri r) 



(45) 



and the function K(t) is the same memory kernel defined in the classical homogeneous 
problem, eq.(8). 
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Given the expressions (41) and (44), the corresponding evolution law (37) for the re- 
duced density matrix is completely determined and provides the general description of the 
dynamics of a quantum Brownian particle acted upon by a multiplicative stochastic non 
homogeneous potential V(x,t) = —Q(x,t)R(t). 



3.b Master equation and decoherence. 



Starting from the functional integral formulation, we can also write the master equa- 
tion for the reduced density matrix, which provides an alternative description of quantum 
Brownian motion. 

A differential evolution equation for the density matrix p(x, x', t) exists if the effective 
Lagrangian corresponding to the effective action (42) is local in time. This happens in the 
white noise limit, in which the influence phase reads 

$[x,x r \ = M 7o J'" dt(Q(x',t) -Q(x,t))(Q(x,t) + Q(x',t)j 

+ Mj J" dt(Q(x', t)-Q(x,t)) (xd x Q(x, t) + x'd x ,Q(x', t) J 



+ 



.M 7o 



-^P jf dt(Q(x',t)-Q(x,t) 



(46) 



The white noise master equation can now be derived following the same procedure 
applied by Caldeira and Leggett in the linear homogeneous case [II]. Starting from the 
effective action (42) with the nonlinear influence phase (46) we obtain 



dp 
dt 



ih ( 2 



d 2 



2M \ dx 2 dx' 2 



%M 



p --iV(x',t)-V(x,t)\p 



^(<2(z',*)-<20M) 

- ^[Q(x,t)-Q(x',t)\ 



'dQ(x,t) + dQ(x',t)' 



P 



\ dt dt 
'dQ(x,t) d dQ{x',t) d 



dx dx 



dx' dx' 



2M 7o 

ph 2 



Q(x,t)-Q(x f ,t)\ p. 



(47) 



The first two terms on the right hand side of eq. (47) represent the mechanical evolution 
without bath-particle interaction. The fourth one describes dissipation and the last one 
the random fluctuations which destroy the quantum coherence. The third term is present 
only when the bath-particle coupling is time- dependent, and it can completely modify the 
usual decoherence patterns as we will clarify below. 

The Caldeira-Leggett master equation in a homogeneous environment is recovered for 
Q = x, that is for d x Q = 1. 



13 



In the classical limit eq.(47) reduces to the Fokker-Planck dynamics discussed in Sec. 2b. 
This can be best understood by transforming the master equation (47) for the density 
matrix in the evolution equation for the associated Wigner function. The task is achieved 
by introducing the "relative coordinate" Y = x' — x and the "center of mass coordinate" 
X = (x + x')/2 and by Fourier transforming eq.(47) respect to the Y variable. In this 
way, one obtains the evolution equation in phase space for the Wigner function W(x,p, t) 
defined as 

W(x,p,t) = JdYp(x-^,X + ^y^ Y / h . (48) 

It is then easy to prove that, in the limit of small H, the evolution equation for the 
Wigner function W(x,p,t) reduces to the Fokker-Planck equation (31) for the classical 
probability density W(x,p,t), where p is the momentum of the particle. 

The quantum model provides very interesting information concerning the dynamics of 
decoherence and the transition from quantum to classical. We consider in particular the 
coefficient g(x,x',t) of p(x,x',t) in the last term of the master equation (47), 

g(x, x', t) = ^-J^ (q(x, t) - Q(x', t^J , (49) 

which is always positive and represents the space-time dependent decay rate of the off- 
diagonal elements in the density matrix; it becomes very large in the semiclassical limit, 
since it is proportional to hT 2 |i~9| , and it reduces to a quadratic function of the difference 
\x — x'\ in the homogeneous case. In the non homogeneous case this term, as well as the 
dissipative one, loses its invariance properties under reflection and translation, because of 
its dependence on the noise modulating function Q(x,t). 

When the function Q depends explicitly on time, eq. (47) describes a decoherence process 
with multiple nonequilibrium scales due to time-modulation of the noise. Explicitely, the 
different scales are represented by the decoherence rate g, the damping rate 7, the external 
potential V(x,t), and the modulating function Q(x,t). 

Overall, these nonequilibrium effects are expressed by the third term on the r.h.s. of 
the master equation; formally, it is similar to the terms derived from a potential, that do 
not involve space-derivatives, with the crucial difference that it cannot be written as a 
difference of the type V(x', t) — V(x, t), and so its effect is not equivalent, as in the classical 
case, to adding a suitable external potential. When the physical situation is such that 
the bath-particle coupling is strongly time-dependent, possibly with very fast variations, 
this term yields a fastly-oscillating contribution in the solution and so can reduce or even 
inhibit the decoherence process. 

Generalizations of eq.(47) may be considered in order to take into account other effects 
such as bath-particle interactions nonlinear in the oscillator coordinates or the finite size 
of the bath. 

In the first case the decoherence rate g(x,x',t) does not depend quadratically on the 
difference Q(x, t) — Q(x', t) (or on the difference x — x' in the homogeneous case), because 
such a dependence is a direct consequence of the oscillator Lagrangian being quadratic. In 
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the second case quantum interference is suppressed only as long as the system is confined 
inside the thermal bath, which is assumed to possess a finite scale £, of the order of its 
linear dimensions. 

The variable Y = x — x' is conjugate to the variable p (see the definition of the Wigner 
function) which in the semiclassical limit reduces to the momentum variable, and actually 
represents the quantum uncertainty on the position of the Brownian particle ||20|| . We then 
expect expression (48) to be replaced by a function g(x, x', t) which becomes negligible for 
\x — x'\ ^> £. 

A finite value of (48) for \x — x'\ 3> £, both in the homogeneous and in the non 
homogeneous models, is possible only if the bath is assumed to be infinitely extended. 
This point clarifies some questions recently raised [21] concerning the decoherence process. 

In this paper we generalized the homogeneous model of classical and quantum dissipa- 
tion in order to include general bath-particle interactions nonlinear in the particle coordi- 
nates. The linear-dissipation theory emerges as a particular case of the model. Moreover, 
we have shown that the nonlinear bath-particle coupling leads to the appearance of non 
homogeneous random and dissipative forces and of a further deterministic time-dependent 
force, in the Langevin equation at the classical level, and in the master equation for the 
quantum case. We illustrated how this yields new interesting consequences both in the 
classical and in the quantum regime. 
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